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Abstract 

In this paper, we study the inversion formula for recovering a function from its 
windowed Fourier transform. We give a rigorous proof for an inversion formula which is 
known in engineering. We show that the integral involved in the formula is convergent 
almost everywhere on R as well as in L p for all 1 < p < oo if the function to be 
reconstructed is. 
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1 Introduction and the Main Result 

The Fourier transform is a very useful mathematical tool, which has been widely used in 
characterization of function spaces as well as in signal and image processing [6J [11] . For a 
function / € L 1 (IR), the Fourier transform of / is defined by 



/>) = / f(x)e~ 



To study local properties of functions (signals), the windowed Fourier transform, also 
known as short-time Fourier transform, is introduced. 

Given a window function g(x), the windowed Fourier transform of a function / with 
respect to g is defined by 



(F g f)(t,u>)= / f(x)g(x-t)e-^dx. 

It is easy to see that F g f is well defined if / G L P (M) and g € L p ' (M), where p, p' > 1 and 
l/p + l/p' = l. 

Continuous and discrete windowed Fourier transforms have been discussed extensively 
in the literature since they are widely used in communication theory, quantum mechanics, 
and many other fields. We refer to IH El IS] for an introduction to the windowed 
Fourier transform. 



*This work was supported partially by the National Natural Science Foundation of China( 10971 105 and 
10990012) and the Natural Science Foundation of Tianjin (09JCYBJC01000). 



1 



Finding a computationally efficient algorithm for the inversion of windowed Fourier 
transforms is a fundamental topic in both theory and applications. The classical method 
to recover / from its windowed Fourier transform is to use the following inversion formula, 

f(- x ) = YTp [f(F g f)(t,u;)g(x-ty™dtdw, (1.1) 

^* I \9 II 2 " " 

R 2 

where we assume that g G L 2 (R). It can be shown that the convergence is in L (R) as 
well as in many other spaces if the function to be reconstructed is and g satisfies some 
further conditions [Tj- 

Since a double integral is involved in (II. ip . it is obviously very complicated. An 
alternate method is to use the filter-bank summation [I], 

f(x) = —L= / {F g f){x^)e^<ko, (1.2) 

where we assume that g(0) ^ 0. Note that (| 1 . 2 j) was presented in [1] in a discrete version 
for compactly supported window functions and the authors stated that their results may 
be equally well stated in a continuous time-domain setting. 

Although (ll.2p is well known in engineering, the convergence of the integral is not well 
stated in literature. In this paper, we show that the integral in (|1.2p is convergent in L P (R) 
for all 1 < p < oo if the function / is. Moreover, by applying the Carleson-Hunt theorem, 
we also show that the convergence is almost everywhere on R. 

Before stating our result, we introduce some definitions. Throughout this paper, xq is 
a fixed real number. For any A±,A 2 > 0, define 

(T AlA J)(x) = f ' {F g f){x - x ,oj)e ix »cLo. (1.3) 
J-Ax 

Our main result is the following. 

Theorem 1.1 Suppose that g is continuous and that g,g £ L 1 (R). Then for any f £ 
L P (R) ; 1 < p < oo, we have 

lim \\T AlA J - 27r^M/|| p = (1.4) 

A\,A 2 — >oo 

and 

lim (T A f)(x) = 2irg{x Q )f{x), a.e., (1.5) 

A— >OD 

where we use the shortcut T A f = T A<A f. 

Remark 1.2 The reconstruction formula is stable in the sense that for any /, / € 

L p (R), 

\\T Al Mf ~ f)\\p < 2 ^ll5lli||/ - /IU VA lt A 2 > 0, 
where C p is a constant depending only on p. For details, see the proof of Theorem \l.l[ 

In Section 2, we give the proof of Theorer dl.il which is based on the famous Carleson- 
Hunt theorem [21 [9] for Fourier series and the extension to Fourier integrals by Kenig and 
Tomas [TO]. 
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2 Proof of the Main Result 

In this section, we give the proof of the main result. 

We begin with a simple lemma on the Fourier transform, for which we omit the proof. 

Lemma 2.1 For any f G L 2 (R) with f G L 1 (M), we have 

/(*) = ^-f f(u>)e ix »dw, a.e. 



2vr _ 

We also need the following formula on the windowed Fourier transform. 
Proposition 2.2 ([8, Lemma 3.1.1]) For any f,g G L 2 (IR), we have 

{F g f){x,uj) = ±-(F- g f){u,-x)e-^. 

Next, we show that for / G L 2 (R) with / G T Al A2 f is convergent in L°° norm. 

Lemma 2.3 Suppose that g is continuous and that g,g G L^R). Then for any f G L 2 
with f G L 1 (M) J we have 



lim \\T Ai ,aJ ~ 27rg(x )f\\ oo = 0. (2.1) 

Proof. For any / G L 2 (R), we see from Proposition 12.21 that 

(T AuA J)(x) = \ M (F g f)(x-x ,u>)e ix »cLj 
J-A x 

= ±- [ M (F g f)(u,x -x)e-^- x ^e^ du 
i- / 2 e ix °»doo [ f(yMy-")e- w{x °- x) dy 

-Ai JR 



2tt 



i- / f(vy yx dy \ M g{y-uo)e-^y-^du: 

-Ai 

cy+Ai 



2tt 
1 

2^ 



Jj/-A 2 

' f W)e-™°»<kj f U+M f(y)e^ x dy, (2.2) 



2tt 

where we use Fubini's theorem twice. By Lemma |2.1| for almost every x, 



{T Al ,A 2 f)(x) - 2irg(x )f(x) 

■w+A 2 



g{uj)e- lxoul doj 



{h[l A 2 hyy^dy-fix)^ 



[ g^e^du [ f{y)e^ x dy. 

or y ~>us-\- A<2 



Hence 

\\T Ai ,aJ - 2vr 5 (xo)/|| 00 < ^ / \g(u)\du> [ \f(y)\dy. 



or B>o)+A 2 



By the dominated convergence theorem, we get 



lim \\T AlA J - 2vr 5 (x )/||oo = 0. 

This completes the proof. □ 
In the followings we prove the convergence in L P (M). First, we show that T Al A2 1S wen 
defined on U>(R). 

Lemma 2.4 Suppose that g is continuous and that g,g € L 1 (R). For any A\,A2 > 0, let 

( sin Ai(y — x) + sin ^b(y — x) 



K AuA2 (x,y) = g(y-x + x ) 



2 sin (^-*0fa-*) S in 



(y - s) 

T/ien we /iai>e 



(2.3) 



(TA u A 2 f)(x) = / f(y)K AuA2 (x,y)dy, V/ G 

./TO 



Proof. Since £ L 1 (1R), we have 5 G for all 1 < p < 00. Hence i 7 ^/ is well 

defined for any / € L P (]R). We have 

(T AuA2 f)(x) = (F g f)(x-x ,u)e i ^du; 
J -A x 
<-A 2 



duj / f(y)g(y-x + x )e-^e™"dy 

-Ai JR 

dy ( M f(y)g(y-x + x )e-^- x ^du J 
J-A x 

f(y)K Al ,A2( x ^y)dy, 

where Fubini's theorem is used. This completes the proof. □ 
The pointwise convergence of Fourier series is a deep result in harmonic analysis. 
Carleson proved that the Fourier series of a function in L 2 [—ir,ir} is convergent almost 
everywhere [2]. Hunt [9] extended this result to L p [— ir,ir] for 1 < p < 00. And Kenig 
and Tomas |10| proved the pointwise convergence of Fourier integral on L P (R). For our 
purpose, we cite the Carleson-Hunt theorem in the following form. 

Proposition 2.5 For A > and 1 < p < 00, define 

(S A f)(x) = [ f(y) Sm f iX ~, y) dy, f € L P (R). (2.4) 
Jr n{x - y) 

Then Sa is a bounded linear operator on L P (M) and there exists some constant C p such 
that 

"sup|(5 A /)(x)| <C p \\f\\ p . 

A>0 P 
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The Fourier multiplier is a useful tool in the study of Fourier transform. The following 
result on the Fourier multiplier is useful in studying the convergence of Ta 1 ,a 2 ■ 

Proposition 2.6 ([6, Corollary 3.8]) Suppose that h is a function of bounded variation 
on R and that (TfJ = h ■ f for f G L 2 (M). Then T can be extended to an operator on 
L P (R), 1 < p < oo and 

\\Tf\\ p <C p V h \\f\\ p , V/gI7(R), 

where Vh is the total variation of h on R and C p is a constant depending only on p. 

The following lemma shows that TA ll A 2 f converges to / in L P (M) whenever / is in 
C*(R), the space of all continuous differentiable functions which are compactly supported. 

Lemma 2.7 For any f G C*(K), we have 



lim \\T AuA2 f-2irg(x )f\\ p = 0, 1< p < oo. (2.5) 

Ai,A2—>oo 

Proof. Fix some / G C*(R). Suppose that supp/ C [ — S~2, for some constant SI > 0. 
Since /, /' G L 2 (R), we have / G L X (R). By Lemma E31 



lim \\T Ai ,aJ ~ 27rg(x Q )f\\ 00 = 0. (2.6) 

Ax,j4a->oo 



Next we assume that 1 < p < oo. By (|2.6p . we have 

lim \\(T Al ,A 2 f ~ ^g(x )f) • X[-2n,an]llp = °- ( 2 - 7 ) 



Al,A2~ >OQ 

On the other hand, put 



K (x,y) = ^ 9 ^°° x^y. 
\x - y\ 



By Minkovski's inequality, we have 



v \ i/p 

|x|>20 ^|y|<n 

(\ i/p 
/ |K(x,y)Rx] dy 
J\x\>2VL ) 

< 4NU f \M\(f T^-z) ' dy 

J\y\<U V J\x\>2U F _ V\ V ) 



M PII/IU 

f p (2Q) 1 1 1 _/ up, 
where M p is a constant and 1/p + l/p' = 1. Note that 



< M p (2fi) 1 / p ' 



\(T AuA J)(x)\< [ \K AuA2 (x,y)f(y)\dy < [ K(x,y)\f(y)\dy. 
J\y\<n J\v\<si 
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By the dominated convergence theorem, we have 

\\( T *i>Aaf ~ 27T 9(xo)f) ■ XR\[-2n,20] Hp = °- ( 2 - 8 ) 

Now the conclusion follows by combining (12. 7p and (12. 8ft . □ 
We are now ready to give the proof of the main result. 
Proof of Theorem II. 1L First, we prove the convergence in L P (M). 
For any / G L 2 (R), by (pT2j) . we have 

(T Al>A2 f)(x) = — / /(y)e^dy / ^K^ ^. (2.9) 
27r JR Jy-A 2 

Hence 

(T AuA jy(y) = h Al ,A 2 (y)f(y), 

where 

Obviously, /iAi,A 2 is of bounded variation on R and Vft, A 4a < 2||g||i. 

By Lemma 12.41 and Proposition 12.61 T Al)A2 is a bounded linear operator on L P (M) and 

\\T Al , A J\\ <2C p \\g\\ 1 \\f\\ p , V/Gl7(R). (2.10) 

Fix some / G L P (R). For any e > 0, there is some / G C£ (K) such that ||/ - /|| p < e. By 
Lemma 12.7} we can find some Aq > such that for any A\,Ai > Aq, 

\\T AltA J-2n'ifao)f\\ p <e. 

Consequently, 

\\T AuA J - 27r^o)f\\ p < \\T AuA . 2 (f - f)\\ p + \\T AuA J-2ng\%)f\\ p 

+ 2k\ 9 (x )\ ■ \\f - f\\ p 
< (2C p \\g\\ 1 + 2TT\g{x )\ + l)e. 

Hence 

lim \\T AlA J-2ngTx7)f\\ p = 0, V/ G L P (U). 

Next we consider the pointwise convergence. For A > 0, let S A be defined by (|2.4p . 
Then Sa is a bounded linear operator on L P (M) and 

(5 , A/T=/-x h A,A]> /e^ 2 W. 

For / G L 2 f|£ p W, define 

(f A f)(x) = -L f gJJje- ix °"(M^S A M u f)(x)dLJ, 

where the operator is defined by 

(M ul f)(x) = e-^f(x). 
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is dense in L P (R), T A can be extended to a bounded 



Since g G L^R) and L 2 f|- 
linear operator on L P (R). 

On the other hand, for any / G L 2 P L P (R), we see from (|2.2p that 



(Ta/)(x) 



to JR Juj~A 



1 



to J K J -a 

1 



, sH e - M °"(M_ w S A M w /)(x)du; 

Z7T 

Using the density of L 2 f| L P (R) again, we get that T A = f A on L P (R). Hence 



(T A /)(x) = — / g{u)e-™°"{M-uS A MMx)du, V/ G L P (R). 



It follows that 

sup \(T A f)(x)\ < f \g(u)\ ■ Bvp\(M- u S A M u f){x)\du. 
A>0 to J K A>0 

By Minkovski's inequality and Proposition 12.51 we have 



su P |(r A /)(x)i 

A>0 



< — \g(u)\ 
to Jr 



swpKM-uSaMMx) 
A>0 



dw 



1 

to 
1 



swp\(S A M w f)(x)\ 
A>0 



dto 



< — / \g{uS)\ ■ C p \\Muf\\ p dw 



Cp 
2^ 



\9\\i 



V/G 



Fix some / G L P (R). For any e > 0, we can find some / G C^(R) such that 

\\f-f\\ P <e. 

By Lemma \2.7\ we have 

lim \\T A f -2tt gJxTlf \\oo =0. 
A— too 

Note that T A f is continuous on R, thanks to Lemma 12.41 We have 

0. 



lim sup 



(T A f)(x)-2irg(x )f(x) 



Hence 



lim sup \(T A f)(x) - (T A ,f)(x)\ =0, Vx G R. 

A,A'->oo 



(2.11) 



(2.12) 



It follows that 



limsup \ {T A f)(x) - (T A ,f)(x) 

A,A'->ao 



< 



lim sup 

A— >oa 



+ 



< 2 



+ 



sup 

A>0 



(T A (f-f))(x 

(T A ,(f - f))(x) 
(T A (f-f))(x 



lim sup 

A,A'~>oo 



(T A f)(x)-(T A/ f)(x) 



lim sup 

A' -too 



< — ll5lli||/-/ll P (using 

7T 



< 



C n 



7T 



\9\\l ■ £• 



Since e is arbitrary, we have 



limsup|(r A /)(x)-(r^/)(x)| 

A,A'-^oo 



0. 



Hence the limit lim A ^ OQ (T A f)(x) exists almost everywhere. Since T A f tends to 2irg(xo)f 
in L P (M), we have 

a.e. 



lim (T A f)(x) = 27rg(x )f(x), 
A— too 



This completes the proof. 



□ 
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